Transport and Strong-Correlation Phenomena in 
Carbon Nanotube Quantum Dots in a Magnetic Field 
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Transport through carbon nanotube (CNT) quantum dots (QDs) in a magnetic field is discussed. 
The evolution of the system from the ultraviolet to the infrared is analyzed; the strongly correlated 
(SC) states arising in the infrared are investigated. Experimental consequences of the physics are 
presented — the SC states arising at various fillings are shown to be drastically different, with 
distinct signatures in the conductance and, in particular, the noise. Besides CNT QDs, our results 
are also relevant to double QD systems. 
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Since their discovery, carbon nanotubes (CNTs) have 
been the subject of intense activity; [l| in particular, ex- 
periments on transport in CNTs have revealed a wealth 
of exciting phenomena. Indeed, long metallic CNTs have 
been shown to behave as quantum wires; 0, H[ negative 
differential resistance has been observed in semiconduct- 
ing CNTs.Jl] Furthermore, short CNTs have been shown 
to behave as quantum dots (QDs),[3j, |5L |6[ exhibiting 
Coulomb blockade (CB) phenomenology [7( known from 
gated two-dimensional semiconducting structures. 

QDs have spurred a renewed excitement about the 
Kondo effect (KE),as they allow detailed investigations 
of the phenomena. [8( In this regard, CNT QDs are ideal 
for studies of Kondo physics. Indeed, initial experiments 
displayed an SU(2) KE arising from the electron's spin; [6] 
more recently, orbital [p as well as SU(4) KEs have 
been observed, d [HI [n( Furthermore, CNT QDs afford 
the possibility of tuning between a variety of strongly- 
correlated (SC) states with a magnetic field. Hoj 

In this work, we consider transport through CNT QDs, 
focussing on their behavior in a magnetic field. We ana- 
lyze the system's evolution from the ultraviolet (UV) to 
the infrared (IR) fixed points (FPs); we discuss the KEs 
that arise and their consequences. More specifically, we 
consider the KEs arising from a single electron (referred 
to as 1/4- filled) as well as two electrons (referred to as 
1/2- filled) occupying the energy levels of the CNT QD 
closest to the Fermi energy Ep of the leads. While previ- 
ous works detailed the properties of the 1 / 4- filled QD.fl2j 
we show that the KEs arising from the 1 / 4- filled and 1/2- 
filled QDs are drastically different; these differences have 
pronounced observable consequences. 

In what follows, we will be interested in the system's 
low-energy physics; hence, we focus on the energy levels 
of the CNT QD closest to Ep of the leads. In the ab- 
sence of magnetic fields, there are two degenerate energy 
levels .[9l[T3| which we label as a and /3. The Hamiltonian 
we consider is 



#QD = ^ (^-iV y-^]T(n QS -rW (1) 

S 

+ E ( [*i ^(°) + h V>L(0)] d KS + h.c.) , 
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where ifil KS (0) creates an electron (at x—0) with spin-s 
in band-K from lead-i (i=l,2); d\ s creates an electron 
with spin-s in orbital-K (n=a, (3) on the QD; n KS =dJ. 8 d K3 



and N=J2 K s " es ! -^o i s tne optimal number of electrons 
on the QD, which can be controlled by a gate voltage; 
Ec is the charging energy; ti is the tunneling matrix 
element between lead-* and the QD; ho is a magnetic 
field. In this work, we take the {ti} to conserve the orbital 
quantum number (which is relevant to the experiments 
m Refs. U3 and [Il|); [3 as a result, the system has an 
SU(4) symmetry when /i =0.[l5j] ho, which would arise 
from a magnetic field applied parallel to the CNT's axis, 
splits the a and (3 orbitals. Throughout this work, we 
employ units where h=l. 

It should be noted ho would also give rise to a Zee- 
man splitting, but this splitting is considerably smaller 
than the orbital splitting, particularly for larger diam- 
eter CNTs. Indeed, the orbital moment /x or b of a 5nm 
diameter CNT was found to be /x or b— 1.5meV/T[13] i.e. 
/iorb— 26/is- [fJ-B is the Bohr magneton.] As we will be 



interested in small fields — /io~0(Tj 
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is given by Eq. [5] — the Zeeman splitting will have very 
small effects. Therefore, in what follows, we focus on the 
orbital splitting. 

We begin our discussion of the properties of CNT QDs 
by considering the current I=(I), where I is the current 
operator 



I = 



^iE[C(oK s -4ii M (o)] (2) 



(e is the electron's charge); in particular, we compute the 
conductance G=dI/dV vs. No (in linear response). We 
are interested in the behavior of G as ho is varied, as well 
as how G evolves (with temperature) from the UV to the 
IR FPs. To understand the IR behavior, G was com- 
puted as per Ref. [H using the logarithmic-descretization 
embedded cluster approximation (LDECA)[l7j] and the 
Friedel sum rule[i~8|: to treat the UV regime — T^>T i7 
where r i =27r ( o i 2 with p being the electrons' density of 
states in the leads — we employed a master equation 
approach. [l9j] 

Fig. QJ, shows G vs. iVo in the UV regime for several 
values of ho. Letting r =2rir 2 /(ri+r 2 ), 



G = e 2 T 
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FIG. 1: G—dl/dV vs. iVo in linear response for several values 
of ho: (a) T=Q.1E C and (b) T=0. 



where PncN^ is the probability for the QD to be in 
a state with N a (Np) electrons in the a ((3) orbital, 
ENaNp is the energy of the state with M^ a N fj being the 
number of these states, and the {N a , Np, N' a , N'^} sat- 
isfy (N a +N )-(N^+N'^)=l. In Fig. [TJi, we observe the 
well-known CB peaks for No=N+l/2 (N is an integer) 
and valleys for No—N. When ho=0, the system has an 
SU(4) symmetry; the two middle peaks have more spec- 
tral weight e.g. the peak at N =3/2 (due to fluctuations 
between states with N=l and N—2) has more spectral 
weight than the peak at N =l/2 (due to fluctuations be- 
tween states with iV=0 and N=l). From the above ex- 
pression for G, this occurs because there are more states 
with N=2 than iV=l or N=0. When h ^0, the SU(4) 
symmetry is reduced to SU(2); as a result, the peaks are 
split and the spectral weight becomes evenly distributed. 

Fig. Qp shows G/G vs. N at T=0, where 
Go = (e 2 /7r)4rir 2 /(ri+r 2 ) 2 . Rather than four peaks, we 
see three distinct plateaus when ho—0 — G/Go—1 for 
the plateaus centered about Nq=1 and No=3; G/Go=2 
for the plateau centered about Nq—2. Furthermore, ho 
has interesting effects on G — whereas ho mainly splits 
the peaks in the UV regime (Fig. [T^), ho has more dras- 
tic effects in the IR. Indeed, the plateau centered about 
Nq=2 is suppressed by ho; the plateaus centered about 
Nq=1 and Nq=3, on the other hand, are unaffected. As 
discussed below, the behavior at T=0 occurs because SC 
states between the QD and leads are formed; ho has dras- 
tic effects on the SC states. 

We now address the physics behind Fig. [1] — we in- 
vestigate the SC states which arise in the IR, as well as 
how they evolved from the UV FP. To this end, we ex- 
amine the QD's spectral function (SF), Ad(u>). Fig. [5] 
shows A d (oj) (at T=0) obtained via the LDECA. For 
comparison, results for Ad(ui) at the UV FP — obtained 
by formally setting {U}=0 — are shown in the insets. 
Fig. [2k shows Ad{uj) at the No— 1/2 CB peak. Here we 
see a broad peak near lo—0 i.e. near Ep ; its features 
do not change much with ho- From the insets, we see 
there was a redistribution of spectral weight, with much 
of the peak's weight in the IR having been transferred 
from higher energies. 

Figs, dp and [2b show Ad(u>) in the CB valleys. A key 
feature is the narrow resonance which appears at or near 
Ep — the Kondo resonance (KR). This resonance is a 
consequence of the SC state formed between the QD and 




FIG. 2: QD's spectral function A d (uj). (a) A d {u>) for jV =l/2. 
(b) Ad(u>) for No=l. (c) Ad{ui) for No=2. Insets: Comparison 
of Ad(ui) at the UV (gray dotted lines) and IR (solid black 
lines) fixed points. 



leads due to the KE; its width represents the dynam- 
ically generated scale characteristic of the SC state - 
the Kondo temperature, Tk- As discussed below, the 
position and width of the KR are characteristic of the 
particular Kondo fixed point (KFP). 

Fig. Wp shows Ad(u)) in the Nq—1 valley i.e. the 1/4- 
filled QD. For /io=0, Ad(ui) exhibits a KR near Ep; for 
/iot^O, the resonance moves toward Ep and its width nar- 
rows. As mentioned above, when ho=0 the system has an 
SU(4) symmetry; ho^=0 reduces this symmetry to SU(2). 
For the 1/4-filled QD, the system flows to the SU(4) KFP 
when ho=0, while /io^0 drives the system to the SU(2) 
KFP; the KR is near (at) Ep at the SU(4) (SU(2)) KFP 

with T£ U(2) <T£ U(4) . The UV and IR behaviors of A d (u>) 
are compared in the insets — the KR is indeed an IR 
property, with its spectral weight taken from the higher 
energy UV peaks; interestingly, ho does not change the 
qualitative features of A d (ui) at either the UV or IR FPs. 

Fig.H: shows A d (oj) for N =2 i.e. the 1/2-filled QD - 
its behavior is drastically different from the SFs arising 
for both N =l/2 and N =l. For h =0, A d (uj) exhibits a 
narrow KR at Ep; for ho^O, the resonance splits and is 
suppressed. Hence, contrary to the 1/4-fillcd QD where 
h drives the system from the SU(4) KFP to the SU(2) 
KFP, h destroys the KE for the 1/2-filled QD. The UV 
and IR behaviors of Ad(to) are compared in the insets — 
we see the KR suppressed as ho increases; as this occurs, 
the peaks at lo=±Ec /2 regain spectral weight. 

Having discussed the QD's SF in the various regimes, 
we now discuss (further) consequences of the SF's fea- 
tures in the CB valleys i.e. for No—N . To facilitate the 
analysis, we integrate out charge fluctuations on the QD; 
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we arrive at the Coqblin-Schrieffer Hamiltonian[18j 

H Q n = ~ (VL /„.) (/«v^„v)-y /kVL'4< 5 ( 3 ) 



where ^.=Mi,JO) + * 2 ?A 2ks (0)]A with t=^f+t|, 
J = (4i 2 /£ c ) [(iV - N - 1/2)- 1 - (AT - N + 1/2)- 1 ] , 
and the fermion operators satisfy the constraint 
fl, s f Ka =N with N being the number of particles on the 
QD. (While discussing the physics of the CB valleys, we 
write the QD's fermion operators as {/ KS }; also, Einstein 
summation convention is utilized). To treat Eq. [3l we 
consider a path integral representation of the partition 
function - we enforce the constraint P KS f KS =N with a 
Lagrange multiplier field A; we decouple the Kondo in- 
teraction using a Hubbard-Stratonovich field x-[HI We 
arrive at an effective Hamiltonian 

tfcff = -y fl< K >f K , s + A (flf KS - N) (4) 

+ j\x\ 2 + X^fl*P K s + X^lf KS ■ 

We begin by considering the physics at higher energies, 
focussing on the flow from the UV to the IR FPs. To do 
so, we treat the Bose fields x an d A in Eq.[4]in mean-field 
theory (MFT). Treating A in MFT amounts to treating 
the constraint P KS f K8 =N on average: {fLfK S )= N - To 
describe the physics near the UV FP, we take (x)=0; the 
physics of the KE is contained in the effective action for 
X, obtained by integrating out the /-fermions and leads. 
To one-loop order, the propagator of the x held J{iui m ) 
is given by the diagram in Fig. [3] (with iv m being a boson 
Matsubara frequency) . Physically, J(iu> m ) is the running 
Kondo coupling. [20l| 

Using our result for J(iw m ), the current I={I) was 
computed as per Ref . fl6l; results for G/Go vs. T/T^ 3 ^ 
are shown in Fig. [H where 



T, 



SU(4) 



K 



Dexp (-1/poJ) 



(5) 



with D being half the leads' bandwidth. [As before, 
G =(e 2 /7r)4rir 2 /(ri+r 2 ) 2 .] Fig. U* shows results for 
the 1/4-filled QD. To begin with, we see that G grows 
logarithmically as T is reduced — this is a conse- 
quence of the logarithmic growth of the running Kondo 
coupling. [l8| Furthermore, G/Go always grows to 0(1) 
i.e. the system always flows to strong coupling. This is 
because the 1/4-filled QD exhibits a KE, irrespective of 
the value of ho. However, G grows more slowly for larger 
h — the system flows to the SU(4) (SU(2)) KFP for 
smaller (larger) ho; the slower growth of G for larger ho 

(See Fig. [2b.) 



occurs because J^ u ' 2 ' <J'|t u ( 4 ) 
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FIG. 3: x field propagator J(iw m ) — the solid (dashed) line 
denotes the leads' (/-fermions') Green's function. 
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FIG. 4: G/Go vs. T/rf U(4) in linear response in the CB 
valleys: (a) 1/4-filled QD and (b) 1/2-filled QD. 



Fig. gb shows G/Go vs. T/T^. U(4) for the 1/2-filled 
QD; the results are drastically different from the 1/4- 
filled QD. As discussed above, whereas ho drives the 
system from the SU(4) to the SU(2) KFP for the 1/4- 
filled QD, h destroys the KE for the 1/2-filled QD (see 
Fig- lit); G even becomes non-monotonic. Such behavior 
has been observed in magnetic alloys, where the occur- 
rence of a spin glass phase freezes spin-flip processes and, 
hence, suppresses the KE.[21| Here, larger values of ho 
freeze both spin and orbital processes. More precisely, 
ho cuts off the growth of the running Kondo coupling — 
for the 1/2-filled QD J(T)=J(iu m =0) is given by 



p J(T) = J < In 



2ttT ^ 



SU(4) 



Re 



UirT J 



where ip( z ) is t ne digamma function: |22| as -0(z)~ln(z) 
for |z|S>l, the growth of J{T) is suppressed for ho suffi- 
ciently larger than TjjL U . 

Having discussed the flow (in the CB valleys) from 
the UV to the KFPs in the IR, we now discuss further 
the physics of the SC KFPs. As before, we treat the 
Bose fields in Eq. 0]in MFT. Now, however, to describe 
the physics of the SC KFPs, we take (xHxo(^0).[H| 
Hence, A and xo are determined via (fl s f KS )—N and 
Xo+2J(^L/ KS >=0- With ( X )^0, the /-fermions SF is 



2r 



(c-e0 2 + r 2 



(e 1/ 2=A±ft /2) J where T=T K when T=0.[3 

Fig. 5] shows G/Go vs. ho/T^ 5 ^ (computed as per 
Ref.Ei) at T=0. For the 1/4-filled QD, G/G =l regard- 
less of the value of ho- This occurs because there is always 
a KE, T^O — for small ho, one is in the SU(4) Kondo 
regime; for larger ho, one crosses over to the SU(2) Kondo 
regime. For the 1/2-filled QD, on the other hand, we see 
that G depends on the magnitude of ho — G/Go =2 for 
ho=0 and decreases as ho is increased. [Within MFT, 
G^0 for /i =2T^ U(4) .] This is because the SU(4) KE is 
destroyed and, consequently, the KR in the QD's SF is 
suppressed for h sufficiently large. (See Fig.[2fc.) 

Also shown in Fig.[S]are results for the noise (which has 
been shown to be a powerful probe of Kondo physics [23l 
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give G/Gq=1 for the 1/4- filled QD, differences between 
the two can drastically be seen in F — F decreases as 
ho increases i.e. as we move from the SU(4) to the SU(2) 
KFP. This is seen most dramatically when t\=ti where 
F^O as ho increases, but the qualitative behavior of F 
is robust. [See the results for ii=2<2-] Physically, this 
arises because the KR in the QD's SF is at Ep at the 
SU(2) KFP, while it is away from E F at the SU(4) KFP. 
(See Fig.[2j>.) For the 1/2-filled QD, F increases as h Q 
increases, approaching unity as G^O; for ti=t 2 , F^O as 
ho^O. For ho—0, the system is in a SC state, with a KR 
at Ep; as a result F=0 when £i=i2- As ho is increased, 
the SC state is destroyed and the system is driven to 
the weak-coupling regime; hence, F—*l i.e. F becomes 
Poissonian. 25] 



FIG. 5: G/Go vs. h / T K and F vs - in linear 

response at T=0. 



|24| |) at T=0. More specifically, we computed the zero- 
frequency noise 



S(eV) = / dt (i(t)I) - (I) 



(6) 



and, subsequently, the Fano factor F=S/2eI in linear 
response. While the conductance probes the spectral 
weight of the KR, the noise gives information about its 
position. Indeed, while both the SU(4) and SU(2) KEs 



To summarize, we considered the behavior of CNT 
QDs in a magnetic field. We analyzed the evolution of 
the system from the UV to the IR FPs. We discussed the 
KEs that occur and their experimental consequences. In 
particular, the KEs arising for the 1/4-filled and 1/2- 
filled QDs were shown to be drastically different, with 
distinct signatures in the system's transport; we are op- 
timistic our results, particularly for the noise. [26j can be 
observed experimentally. Besides CNT QDs, our results 
are relevant to double QDs and, more generally, to QDs 
with two-fold orbital degeneracy. 
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